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^SJ ' Abstract 

. We investigate the phase space of a quintessence theory governed by a generahsed version 

I of the DBI action, using a combination of numeric and analytic methods. The additional 

degrees of freedom lead to a vastly richer phase space structure, where the field covers the 
full equation of state parameter space; —1 <lo< 1. We find many non-trivial solution curves 
to the equations of motion which indicate that DBI quintessence is an interesting candidate 
Q ' for a viable k-essence model. 
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1 Introduction 



> 

\ The dark energy problem continues to be a sticking point for theoretical physicists. The simplest 

■ solution to this problem is to postulate the existence of a vacuum energy or cosmological constant 
. which agrees with all the current observational bounds [UElElll]. However we are then left with 
\ a secondary problem, namely explaining why the vacuum energy is tuned to such a small value 

without some obvious symmetry to protect it. For many years we have hoped that UV complete 
theories of gravity would shed light on this issue, which is in effect an extremely embarrassing 
IR problem from this perspective. However despite much effort, neither string theory nor loop 
^ ' quantum gravity has shed any compelling light on this issue - although there have been many 

■ interesting proposals. 

An alternative approach is to assume that the cosmological constant is exactly zero, since 
supersymmetry can then be invoked as the regulating symmetry in this case. However one then 
has to account for the fact that low-energy supersymmetry must be broken and an alternative 
explanation for the current expansion must be sought. One way to deal with the latter problem is 
to assume that the dark energy phase is driven by a dynamical field, implying that the equation 
of state is an explicit function of time [71 [8]. Currently this cannot be ruled out by our best 
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observations and therefore remains a possible solution to the dark energy problem. However 
one cannot just consider ad hoc scalar fields coupling to gravity, since the low energy theory 
will still be sensitive to high energy physics. In particular we must ensure that any additional 
scalars are neutral under all the standard model symmetries, and that they do not introduce 
additional fifth forces. Therefore one must search for viable models of dynamical dark energy 
within UV sensitive theories. 

Phenomenological models of our universe have proven difficult to construct within string 
theory, due to technical difficulties arising from moduli stabilization, whereby we assume that 
the extra dimensions of the theory are compactified on manifolds with SU{3) x SU{3) structure 
(in the type IIB case) [5], and orientifolded to preserve the minimal amount of supersymmetry 
in four-dimensions. As a result, embedding realistic cosmology into string theory has proven 
difficult. One area which has been well explored in recent years, is inflation driven by the open 
string sector through dynamical Dp-hianes. This is the so-called DBI inflation [111 [T5] - which 
lies in a special class of K-inflation models. It was originally thought that such models yielded 
large levels of non-Gaussian perturbations which could be used as a falsifiable signature of string 
theory [12]. However subsequent work has shown that this is may not be the case, and that 
the simplest DBI models are essentially indistinguishable from standard field theoretic slow roll 
models [22l EH [26] El- The problem is that the WMAP 5 year data set [2] imposes very tight 
constraints on the allowed tuning of the free parameters in the theory. We are then left with the 
choice of either having large non-Gaussianities but with vanishing tensors, or assume that the 
tensor spectrum will be visible - in which case there is no non-Gaussian signature. The models 
are only falsifiable once these conditions are relaxed. One can get around these conditions by 
considering more complicated models such as multiple fields [16] , multiple branes [181 US US] ; 
wrapped branes [27] or monodromies [2^ - but even here there are still problems with fine 
tuning, backreaction and the apparent breakdown of perturbation theory in the inflationary 
regime |29j . 

In models of dynamical dark energy, on the other hand [3 [SI [U] , the WMAP constraints can 
be relaxed and therefore DBI models may still have some use as an explanation for a dynamical 
equation of state. Moreover this fits in nicely with several intuitive ideas from string theory. 
Namely that inflation can still occur, albeit only through the closed string sector - where one (or 
more) of the geometric pseudo-moduli are actually responsible for the initial inflationary epoch 
(see [6] for the phenomenologically most viable proposals). After inflation the universe lives 
on branes that wrap various cycles within the compact space and are extended along the large 
Minkowski directions. In this sense we see that a GUT or Electro- Weak (EW) phase transition 
can manifest through a geometric fashion - namely the Higgsing of branes in the compact bulk 
space. This suggests that dark energy may well be a dynamical process, and moreover in the 
light of these open string constructions, retains a sense of being geometric in nature. 

With this in mind, various authors have begun to explore the phase space of DBI-driven dark 
energy [13^ 114]. The initial works have dealt with the dynamics of a solitary Z)3-brane moving 
through a particular warped compactification of type IIB. In this note we wish to generalise this 
further to a more phenomenological class of models that include multiple and partially wrapped 
branes. We believe that this may be a more generic situation to consider, since typically one 

^Note however, that the models proposed in [211 126j evade such problems. 
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should expect branes of varying degrees to be wrapped on non-trivial cycles of the compact 
space. Our work is a first step into considerations of a more general set-up for quintessence in 
IIB (open)-string theory, and we hope will be a valuable starting point for further endeavour. 



2 Dynamics of the effective theory 



To begin let us assume that the the universe at such late times can be adequately described by 
a flat FRW metric and the matter sector consists of a dynamical scalar field and a perfect fluid, 
which are both separately conserved. The usual cosmological equations of motion are therefore 
independent of any particular model and can be written as 

Pi = -3HiPi + p,), 

where i runs over the contributing components. The equation of state is given by uji = Pi/pi, 
however if to of the fluid component is assumed to be constant then we can integrate the appro- 
priate conservation equation exactly to obtain 

p oc a-=^(i+") , (2.2) 
where the scale factor varies as a function of time such that a{t) ~ t2/(3[i+aj])_ 



The model dependence arises in the parameterisation of the scalar fleld sector. In our case 
we are assuming that the dark energy is driven by open string modes, which at low energies are 
described by fluctuations of a Dp-brane whose dynamics are governed by the Dirac-Born-Infeld 
action (DBI) - which is a generalisation of non-linear electrodynamics [TS] . Typically one 
assumes that the standard model is localised on an intersecting brane stack, in one of the many 
warped throats that are attached to the internal space. For consistency reasons in the simplest 
cases, these are taken to be either D3 or Z)7-branes. In this note we will consider a bottom up 
approach therefore we shall not worry too much about the geometric deformations of the compact 
space, nor any constraints imposed by Orientifold Op-planes - aside from those that ensure that 
all tadpoles are consistently canceled so that we can trust the low energy supergravity theory. 



The action we consider is a generalised form of the DBI one coupled to Einstein-Hilbert 
gravity, which can be embedded into this background and takes the following generalised form0 



S = -J d^xa^t) (TicP)W{cP)^l-^^-T{<p) + VicP) \ + Sm , (2.3) 

where T{<j)) is the warped tension of the brane and Sm is the action for matter localised in 
the Standard Model (SM) sector. Thus our assumption here is that our dynamical open string 



"^We refer the more interested readers to [23] for more details on the precise structure and origin of this action. 
The important thing to note is that ((> is a, matrix valued field. For recent work in a related direction see [24j . 
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sector is coupled only gravitationally to the SM sector and so we do not have to worry about 
additional forces or particle production. There are two potential terms for the scalar field which 
are denoted by W{(j)) and V{(l))- The first of these terms can arise in different places within the 
theory. Firstly if the brane is actually a non-BPS one [lUj . then the scalar field mode is actually 
tachyonic and the potential is therefore of the usual runaway form. If there are N multiple 
coincident branes, then the world- volume field theory is a U (N) non-Abelian gauge theory and 
the potential term is simply a reflection of the additional degrees of freedom [30]. Through the 
dielectric effect, one can also see that this configuration is related to a D5-brane wrapping a 
two-cycle within the compact space and carrying a non-zero magnetic flux along this cycle. Both 
of these configurations lead to an additional potential multiplying the usual DBI kinetic term. 

The origin of the V{4>) term is less explicit - but is a sum of terms. One expects open or 
closed string interactions to generate a scalar potential V{(j)), however the precise form of such an 
interaction depends upon many factors such as the number of additional branes and geometric 
moduli, the number of non-trivial cycles in the compact space, and the choice of embedding for 
branes on these cycles. Typically one can only compute this in special cases in the full string 
theory. There are also additional terms coming from coupling of the brane to any background 
RR form fields. The action above is assumed to be that of Z)3-brane(s) filling the space-time 
directions, which naturally couple to the field C^^^ through the Chern-Simons part of the action. 
However for wrapped L>5-branes there is also the possibility of a coupling C^^^ A F, where F is 
the magnetic field through the two-cycle. For example in the warped deformed conifold one can 
see that dC^^^ = •kdC^'^'^ and therefore there is an additional term in the DBI action 

S j d'^xa^{t)g-^Ma'T{(j)), (2.4) 

up to a normalisation factor of order one. Terms such as this have been added to the interaction 
potential to define the full scalar potential V^(</'). Recent extensions to standard DBI inflation 
have included the contribution from higher dimensional bulk forms, with the remarkable result 
that they cancel one another up to third order in the action and therefore do not affect the 
leading order perturbations [20]. Extending this work to higher order is therefore extremely 
interesting. 



The corresponding equations for the energy density and pressure of the DBI can then be 
written succinctly as 

P^ = ^b- Wicj))] - VicI)) , = T{cP) [W{cPh - 1] + ^(0) , (2.5) 

where 7 = [1 — (j)"^ /T((f))]~^^'^ is the usual generalisation of the relativistic factor. The subscript 
ij) denotes the scalar field component here. We can also immediately define the equation of state 
parameter for the quintessence field to be 

^ T{cP)[j-W{cP)]-VicP)j 

^ r(0)7[VF(</.)7-l]+n'^)7' 

from which one clearly sees that it is dynamically sensitive and can take a wide range of values. 
For instance we only recover a;,;!, ~ — 1 in the limit that the field is non-relativistic and the entire 
solution is dominated by the V{(t)) terms - which will clearly require large amounts of fine tuning 
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to accomplish. There are clearly several regions of parameter space that are of interest. First let 
us assume that the potential term is zero, either because it is suppressed or there is an unlikely 
cancellation between the contributing terms. The more general case with non-zero V leads to a 
wide variety of complex behaviour. We can therefore identify several limits of interest - focusing 
on the behaviour of W: 



• W{<j)) = 1 - which reduces the action back to the usual DBI case which has = I/7 as 
discussed in [T3]. 

• W{<j)) = 07 - which leads to constant 7 if a is constant, since the two are related via 
7^w<^a = 1 — a + u>^. Moreover this again means that 4> oc t~(-^+'^</')/(^+'^) as in the case 
where W = 1. 

• W{<f)) ^ - as could occur in the case of a tachyonic theory, which mimics a dark energy 
dominates phase with ut^j, = —1. However one must be careful if this is to be representative 
of non-BPS D-brane actions, since the coupling to the form field is non canonical in this 
instance. In fact the coupling term will typically be of the form dcp A C. This means that 
there is no solitary T{(j)) term in the action and therefore the equation of state in this 
instance will vary like — 1 /7^ . 

• W{(j)) 7 - which can occur in the multi-brane/wrapped brane case and yields a;^ ~ 
-1/7^. 



Note that in all cases the equation of state parameter remains bounded between —1 < w,^ < 1. 

One can combine the expressions for the energy-momentum tensor components, and together 
with the continuity equation we obtain the following equation of motion - assuming that the 
scalar field follows a monotonic path 

which is a generalisation of the Klein-Gordon equation for the DBI Lagrangian. The subscript 
of T, W and V denotes derivative with respect to the field value. The other dynamical equation 
of motion for the Hubble parameter can be written as 



H- ' 



2M2 L' 



p{l+uj)+jW{4>)(t>' , (2. 



where we have defined the pressure of the barotropic fluid to he P = ujp and that it is non- 
interacting. We leave the interesting case of interacting pressure for future endeavour. 

Let us consider, as an example solution, the case where there is a scaling solution with W = 1, 
which has been reviewed elsewhere [H]. We will find it convenient to define the quantity 
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in which case we see that 
follows 



This allows us to reconstruct the tension of the brane as 



X = 1 



(2.10) 



where is a dimensionful mass scale, A is a constant and a = 2X/ (1 — X). Using the fact that 
ajfj, = 1/7 we can then see that for X ^ 1 the solution is physically valid only when uj > 2/a 
since we define 7 to be the positive root. Let us now consider the phase-space dynamics of 
the theory in more detail following along the lines of [7]. It is initially convenient to define the 
following new variables 



r(0)VF((/))7 1 

3 HM„ 



7 



1 



3 HMr, 



1^3 



VTMpV^ 



Y3/2 

VTMpT^ 

Y3/2 

W^Mp 
^3/2^5/2 ' 



(2.11) 



in terms of which we can see that 7 = [1 — y^/ (3x^)] ^^"^ and the fluid density parameter can be 
written as 



1 



2 , 2 

J + X 



1 



1 



VF(</>)7 

whilst the equation of state in dimensionless variables will become 

_ 1 / x^[-f-W{^)] - z^VF((/>)7^ \ 
'^'^ ~ 7 \x^[W{^)j - 1] + zm{(l))j ) ■ 



(2.12) 



(2.13) 



As is customary we will now switch to dimensionless derivatives, denoted by a prime, replacing 
time derivatives by derivatives with respect to the e-folding number, M. Therefore we can easily 
determine 



IF 



' 2 



3(1+, 



1 



X 



1 



1 



w{4>)i 



(2.14) 



A useful quantity to calculate is the variation of the kinetic function, which we can write in the 
following manner using the equation of motion 



7 ^ _ _ ^ _ 

7 T W T jWT 



(2.15) 



We can then determine the dynamical equations for the dimensionless fields as derivatives with 
respect to M 



1. .yz'"^ y"^ H' 
7T(w+M2)-^-:T--a;— , 



2 

-32/ ( 1 



2x 



H 
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'IF 



i + — [/ii + H + 

xy 



7X 



3x H3 - y- 



(2.16) 
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and the remaining parametric solutions are 



° W,}' 

X \ 2 2T% J ' 

Note that if the /ij are constants, then the previous three equations form an autonomous set 
and should uniquely specify the dynamics of the quintessence field. We will consider this case as 
the simplest (canonical) example. If we wish to appeal to string theoretic constructions then we 
restrict the parameter space of solutions. It is more interesting to consider the above equations 
in the context of a phenomenological model and see what kind of functions yield the correct 
behaviour. Explicit constructions of string backgrounds are typically difficult and there are only 
a few well known examples that are ritually invoked, however if we take string theory seriously 
then there are undoubtedly other non-trivial backgrounds that are cosmologically interesting but 
not yet constructed. Since an analytic analysis of this generalised system is highly complicated, 
it is convenient to use a combination of analytic and numerical methods to understand the 
dynamics of the system. For a numeric analysis it is necessary to re-write the fiuid equation in 
terms of more useful variables. It turns out that the simplest variables to use are the following 

= (2.18) 

3 



^ _3^ 3Mp£3/V^ 
72 ^ X V 2 



2 1 



ij^^\_3M^_^H^ (2.19) 
^^5/2 J ^JW^ H ' ^ ' 



which are easily derivable from the terms written above. The equations (j2.14p . (|2.16|) . (j2.17p . 
()2.18p and (12.190 together with barotropic fiuid equation: p' = — 3p(AA)(l + w)^ hence form a 
closed ten-dimensional autonomous system if T,W ox V are given as explicit functions of (j) or 
as constants. 



2.1 Case I 



Let us take the canonical string theoretic example arising when the local geometry can be 
approximated by an AdS space. This geometry typically arises in the near horizon limit of 
coincident L'3-branes (or fiux). In this case we see that (at leading order) 

n<^) = ^, y{'p) = '^. w{<t>) = w, (2.20) 

where we have also included an effective (j)^ potential for the system. This means that /^s = 
and we also have a constant /Ui which allows us to write the remaining terms as 

2V2Mp 2x2^1 

W = 72^, /^2 = -77r^, (2.21) 
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and therefore the dynamical equations reduce to 



2x2 
-3y fl 



1 



2x 



xH' 



H 



2x 



zH' 



1 + 



xy 



2x^ 



H 



Q^izx yH' 

~^ IT 



(2.22) 



The simplest way to proceed with the analysis is to consider the final equation above, since this 
splits the solution space neatly into two components. Thus we search for solutions where either 
z = or z = {2x/y^i)H' / H as initial conditions. 



The first sub-set of solutions admits (0, 0, 0) as a (trivial) fixed point, which is a fiuid dom- 
inated solution since = 1 in this instance. Let us remark here that this fixed point solution 
will occur for all the cases we consider, however since this implies a vanishing of the action, 
causality implies that this fixed point must be unstable - i.e. phase space trajectories will fiow 
away from it. By making this field a phantom scalar, one can evade this causal bound and the 
point can become a stable fixed point. This behaviour arises in many places in the literature, 
so we will not discuss it further here. 



There is also a critical point at (1, \/3, 0) which is a kinetic dominated solution. This solution 
actually exists as solutions to the quadratic expression = '^x'^ which corresponds to the limit 
7 ^ oo. In terms of the density parameter, a quick calculation shows that along the general 
curve (parameterised by yo and xq) we find = 1 — Xg. Thus at the trivial fixed point we 
see i7 — > 1, however for xq — > 1 we see that ^2 ^ corresponding to non-relativistic matter, 
i.e. dust. In this instance we also find a{t) ~ t^/s expected from the cosmological evolution 
equations. Again due to the special algebraic properties of the DBI action, we anticipate that 
this solution will also be found for the other cases of interest. 



The second sub-set of solutions are more interesting, as initially one can solve the system 
by slicing the phase space at y = H. One can use the condition on H' to fix z through 
2;2 = 1 — x'^{W — 1)/W. Combining this with the equations of motion gives us the following 
fixed point (taking positive signs of all roots for simplicity) 



y = 0, z = l, (2.23) 

which is valid for all < 1 in order for these points to be real and at finite distance in phase 
space. If we then compute the density of the fluid we flnd = 0, since il^ = 1, which corresponds 
to a purely dust-like solution. Note that this class of solutions does not exist for the simple D3- 
brane analysis as in [U], since it arises from additional degrees of freedom which are neglected 
in these models. The remaining solutions in this sub-set are difficult to find analytically. 

More generally we can see that the above solution corresponds is a special case of the more 
^Note that one cannot do this for x = since the action becomes singular and ill-defined. 
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general Case I behavior, which we paramaterise by 



m 



^('^) = 5^' ^('^) = ~5^' W{4>) = W, (2.24) 
where we can then exphcitly write 

where ^4 is a (real, positive) constant provided that /3 > 0. 

but which simplifies in the limit a = j3 + 2. As before, the solution space splits into two 
disconnected sub-sets, therefore in the first instance where we take slices through 2 = 0, we find 
the following bound 

^^>0, (2.27) 
(a - P) 

which implies that a > [3 and so the brane tension should dominate the dynamics (in the large 
field regime). Let us therefore assume that a, /? are chosen such that this condition is satisfied 
- then we find the solution branch is governed again by the relation = Sx^ as expected - 
which contains the solution (0, 0, 0) as a special case. Moreover this is valid for all values of a, (5 
satisfying the above constraint. The secondary solution branch occurs when we find solutions 
to 

which is generally very complicated. A simple set of solutions do arise when we consider slices 
at 2/ = 0, since the fixed points are localised along the curve 



y = 0, z = ±,/-^, (2.29) 



\l {a- [3){l-Wy ^ ' V(a-/3)' 

which corresponds to a dust-like solution 17 = Va, [3. The reality constraint here demands that 
a > (3 which in turn fixes W <1. However there are also additional solutions where /3 < and 
positive a - provided that W > 1. Explicit realisations of this scenario within a string theory 
context can arise through potentials arising from brane/anti-brane interactions and is therefore 
a non-trivial and interesting solution. 



Figs. [T] and [2] show the numerical solutions in phase space. For W = 1 case, the numerical 
constants are given as Mp = l,m = l,A = l and w = (dust case). Other parameters are 
a = 4, /3 = 2 and A = 2\/2. As expected the (five) fixed points all lie along the curve 7/^ = 3x^. 
We also plot the evolution of each parameter (x,y,z) as a function of the e- folding number in 
Fig. [3] where each of the coordinates tends to its critical value. As expected the phase space 
dynamics are Z2 symmetric about the origin. Note that in the case of y{N) one can keep y 
suppressed for a few e-foldings with enough tuning, before eventually it evolved towards the 
points ib\/3 at late times. The full numerical solution of the case VF < 1 is illustrated in Fig. 
H] where W = 0.95, which uniquely fixes the critical points to be a; = ±^/20,y = 0, z = 1. 
As one can see from the resulting plot, this is an unstable node because the general behaviour 
is divergent. Note that x ^ 00 in this regime effectively solves all the dynamical equations 
trivially. 
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Figure 1: (Case I) 3-D xyz phase space trajectories for T{(f)) = cj)^ / X"^ , V ((j)) = nn?(jP'/2 and 
W{(j)) = W . We have set here, Mp = 1,W = l,m = 1,X = 1 and w = Q (dust case). 
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-3 1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 M 1 1 I r I I 
-15 -10 -5 5 10 15 

X 

Figure 2: (Case I) Phase space trajectories in xy plane. Four attractors (±1, ib\/3, 0) and one 
unstable node (0, 0, 0) can be seen here, z is bounded within (—1,1) range. 
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Figure 3: (Case I) Evolution of x, y, z versus e-folding number setting W = 1. 
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Figure 4: (Case I) Evolution of x,y,z versus e-folding number setting W = 0.95. All solutions 
diverge from the origin 
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2.2 Case II 



Analogous to the first case, let us now consider another branch of solutions where this time the 
tension of the brane is taken to be constant. This dramatically alters the relativistic rolling of 
the scalar field since the 7 factor is no longer warped. Initially let us consider the ansatz 



V{<t>) = ^, T{ct>)=T, ^0) = ^, (2.30) 



which implies that 



^1 = I X4^3 ) ' /^2 = 0, /^3 = ^, (2.31) 



and the corresponding field equations become 

, a'^yz^ xH' 



X 



2x H 



\ bx^y \ Ix'^y J ■jx H 
where we have defined a as the constant pre- factor in the definition of ^1. 



As before we separate the solution space into two - first finding solutions to 2; = and then 
solutions to H' / H = a'yz^ /{2x^). In the first case is is straightforward to see that there are 
the usual fixed point solutions at (0, 0, 0) and (1, Vs, 0) (with their respective partner solutions) 
respectively coming from the usual condition that = Sx^. The secondary branch of solutions 
also admit fixed points when y = 0, however the condition on z is that z = 0, — 4x^. Since we 
want real solutions we are forced to set z = as a secondary constraint. This forces W to diverge 
and therefore in the limit that 2; ^ we find that ^ ±1 which is a unique solution. Again 
the density parameter vanishes identically in this limit as one would expect. The remaining 
solutions are actually extremely difficult to solve analytically as they correspond to high order 
polynomials. As a result we are forced to sketch their behaviour numerically. 



Phenomenologically we see that the ansatz presented above is a special class of the more 
general solution 



r(</.)=r, v{4>) = ^, w{^) = j^, (2.33) 



which has the parameterisation constraints 

(2+/3)/(a-/3) 



a(^) , /.2 = 0, ;.3 = i^7(— (2.34) 



where A, B are both constants. One can see from the dynamical equations that fixed points 
with z = can only occur when the following condition is met 

2(1 + a) -B , , 

' ^ >Q, (2.35) 
a — p 
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which is trivially satisfied for cases where a > (3 (which we assume as an initial constraint). 



More generically we see that, provided a > —2, we recover the usual fixed point equation 
= 3x^. However we need to be careful here because if this condition is satisfied then W 
becomes undefined. Since this is the overall pre-factor multiplying the DBI action, the action is 
undefined in this limit and it should therefore correspond to a point of instability in the phase 
space. In the limit where a = —2, which implies that /? > —2, the fixed point solution now lives 
on the zeros of the polynomial 

3x^i?7-5(2+/3)/2(i+/3) ^ 3^2y - / = , (2.36) 

which can be used to fix x = x{y) or vice-versa. This solution is actually indicative of a more 
general branch of physical solutions where we take /3 > 2(1 + a). The resulting fixed point 
equation (provided a 7^ 2) is trivially calculated to be = as before, but now we see that 
W vanishes identically. In turn this means that the kinetic terms also vanish and the solution 
is dominated solely by the potential interaction. One could imagine a situation such as this 
occurring in the condensation of an open string tachyon mode on a non-BPS brane, where the 
vanishing of W indicates that we are living in the closed string vacuum. For dynamic solutions 
it seems reasonable to consider this particular case as the late time attractor for the solution 

The second sub-set of solutions is again complicated, but again we can analytically under- 
stand the plane at y = 0, which gives us the fixed point solutions 



X = ±\ , y = 0, z = ± 

V a 

Clearly for the solution to be real we require that a, (3 have opposite signs. This satisfies our 
primary constraint, therefore is a physical possibility. Moreover in the limit where we set (3 = —a, 
we find that = which is again the dust solution. Illustrations of numerical solutions for the 
case II are in Figs. [U [6] and [71 Constants are set as Mp = 1, T = l,m = 1, A = 1 and w = 
(dust case). Other parameters are a = 4, /? = 2. From the numerical analysis one sees that 
there are six saddle nodes, only two attractors and one repulsive point which is the origin (0, 0, 0) 
as expected. The dynamical trajectories are particularly interesting due to their apparent lack 
of monotonicity as a function of e-fold number. The z term in particular appears to have a 
large variation in trajectory, diverging in some instances whilst rapidly reaching zero in other 
instances. Conversely the y variable displays very uniform (physical) trajectory behaviour, with 
several curves almost on top of one another at y = and the remainder smoothly driven to the 
(unstable) critical point yc ^ 1.8 in the example given. 
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Figure 5: (Case II) 3-D xyz phase space trajectories for T(0) = T, = m^(/>^/2 and W{(l)) = 
4)^ /X^. We have set here, Mp = 1, T = 1, m = 1, A = 1 and w = (dust case). 




Figure 6: (Case II) Trajectory shce through the yz plane. 
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2.3 Case III 



Let us now consider a new case where only W = W{(j)), with all the other terms being constant. 
We will take W = (j)"/X" for generality - which in turn should impose a constraint on the allowed 
values of a. In this case we see that 

/Z\(a+2)/a 1 

;,i = /X2 = ^^ = aA(-) -T^—TTT- (2.38) 



xJ -)/{2a"l)/" 

where A is a function of the constant parameters A = Mp/A(T/1/)("+^)/^". Because only /Z3 is 
non-zero the resulting dynamical expressions are considerably easy to work with 



x' 



xH' 

y' = -32/(l-^) -3a^z(-+2)/"^{"-2)/a_y|: 
zH' 

z' = (2.39) 

Considering the slice again through z = 0, we see that the solutions split into two types depend- 
ing upon the integer a. We recover the usual = curve only when a > or when a < —2. 
If a = —2 then the corresponding polynomial equation becomes 

2/7*^/2 = 2Ax^ (2.40) 

which is difficult to solve analytically due to the dependence of 7 on both x, y. This expression 
does not admit anything but the trivial solution if we set y to zerclfj. Again we see that there is 
a potential problem here since the potential W goes like l/z^, and is therefore divergent in this 
limit. Solutions to this expression are possible, but complicated. Interestingly there does exist 
a solution curve given by 

22 81 \/3a , . 

y = ax^, Xr = — 7- -TTTT (2.41) 

f c> 2^(9-3o)9/4 ^ ^ 

where the parameter a factor must satisfy < a < 3 for this solution to be physical. Since a 
need not be integer, there are essentially a continuum of curves giving rise to fixed points in this 
theory. 

The secondary branch of solutions again admit fixed point behaviour for y = 0, however 
things are more complicated since the fixed points are now obtained by solving more non-linear 
expressions. There are two cases of immediate interest however. Firstly if we have a = 2 then 
we see that = —1/(2^4) which is only real when A is negative. Since we have chosen our 
parameterisation such that this quantity is positive, this particular branch of solutions is ruled 
out. Interestingly when a = —2 there is a unique fixed point located at 



. = . = 0, . = ±^=L=yi^-l. (Z42) 



'By trivial we mean the point (0, 0, 0) 
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which corresponds to a positive definite equation of state parameter 



^^_ 2T^A{A-l) + V^{T/V-l) 
ATV{T/V - 1){2A - 1) ' 

Note that we must require T > V for this solution to be non-singular, which means (again) 
that the tension term dominates the energetics of the theory. What is also obvious is that 
demanding A = 1/2 leads to a novel fixed point at (±1, 0, 0) regardless of the ratio T/V . Using 
the definition of A this fixes A = 2Mp and therefore W is vanishingly small unless the scalar 
takes is trans-Planckian. This is manifest in a divergence in the equation of state parameter and 
is therefore unphysical. Therefore we must ensure that A <\/2 implying that A > 2Mp. Since 
this is the largest scale in our theory, one again expects this to be unphysical. 

The more general solutions can be found numerically and corr espond to xl = l + zl{T/V -I) 
where Zq are the characteristic solutions to the non-linear equation 

1 + a^z(2+a)/2 ^ ^ z^{T/V - I)) = 0. (2.44) 



In this more general case we can set T = V without the solution diverging, and we therefore 
find the corresponding fixed point solution is thus given by 

/ I \2/(2+a) 

x = ±l, , = 0, -(-^) (2.45) 

which implies that a is negative. Moreover we see that is again zero here for all physical 
values of a, although there is no additional constraint upon the magnitude of A. Now, we see 
numerical solutions in Figs. [SI E]and[TUJ Constants are set as Mp = 1,T = = 1,A = 1 and 
w = (dust case). Other parameters are q = 1 and A = 1. 
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Figure 8: (Case III) 3-D xyz phase space trajectories for T((?!>) = T, V{(j)) = V and W{(l)) = 
(/>"/A". We have set here, Mp = 1,T = V = l,m = 1, X = 1 and w = (dust case), a = 1. 
Green hues approach an attractor 




I [ I I I I I I I I I I I r I 1 ] I I I r il 

X 



Figure 9: (Case III) Phase space trajectories in xy plane. 
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Figure 10: (Case III) Evolution of x, y, z versus e-folding number. 
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2.4 Case IV 



Following on from the previous class of models, we can find solutions where the scalar potential 
is now constant, using the ansatz. 



v = v, T{4>) =\\] , w{4>) 



(2.46) 



where A, 5 are terms of the requisite dimensionality. From this expression we see that is 
identically zero. It will be convenient to define the following function Q = V\°^5l^ which in turn 
can be used in the definitions of the remaining /ij functions 



ni 



aMp (Qx' 



74/2 
3a- 2 



72:^ 



ni 



(2.47) 



-712 



aMr, 



2{a + (3) ' 

and now the dynamical equations simplify to become 



?^2 



1 + 13 
3a - 2 



x' = 


2x2 


2x 


xH' 
H 


y' = 


-3y(l 






z' = 


zH' 






H ' 







xy 



+ 3;U2 



p/3-Q:^2/3+a 




o 2 yH' 
^3 ~ ~jf 



(2.48) 



The resulting analysis is far more complicated than in the previous cases. Let us again start 
with the simplest solution slices at 2 = 0. The expressions for x' and z' readily simplify in this 
instance, however the equation for y' requires us to be more careful. We see that in order for 
the 2;^//2 term to vanish in this limit we require (2 + 3/3)/(a + /3) > 0. The remaining ^2 term 
only vanishes if this condition is tightened to (2 + /3)/(a + /?) > and the term coming from /^s 
only vanishes if (1 + /?)/(a + /3) > 0. If these inequalities are reversed, for example, then these 
terms diverge in the z — > limit. If we restrict ourselves to well-behaved solutions such that 
a,/3 satisfy the above bounds (either by both a, /3 > or by a > 0, /3 < with |/?| > |a|), then 
we obtain the solution curve = ^ usual. If the parameters a, /3 do not satisfy at least 
the minimal bound, then one can only solve these expressions numerically. 



The only other solution branch occurs when H' /H = 0. This is again a complicated solution, 
however things simplify somewhat when we slice through y = 0, but also tune the solution such 
that a = /?, which gives us 



2VQ 



4Q(x2 - 1) 



X 



2520 



(2.49) 
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and therefore the fixed point solution in this instance is given by solutions of the polynomial 




VQ 



) 



l/(2a) (l+8a)/(4a) 



Z 



= 1. 



(2.50) 



This can actually be solved exactly when a = —1, but numerically for more general a. The 
exact case gives us the following solution 



where zq is given by the term written above. This is a highly complicated solution, but one 
sees that in principle there are many fixed points along the plane (xq, 0,2:9) depending on the 
constants A, 5. One also sees that there is a simple solution when x = 1, since this implies that 
Zq = b'^l^fQ or Zq = 0, the latter again giving rise to the point (1,0,0) which corresponds to 
the non-propagating end point of the brane dynamics. 



Q6-2,/Q + 5^ ±Sy^F(Xjj 
2X5^ 



(2.51) 



F{X,6) 



+ 
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X 



Figure 11: (Case IV) 3-D xyz phase space trajectories for T{(\)) = [cj) / X)"- ,V {4>) = V and 
W{(t)) = {(t)/5f. Here, Mp = 1, F = 1, m = 1, A = 1, a = 1, /3 = 1, 5 = 1 and w; = (dust case) 




+: 

6ll I I I I I I I I I 11 I I I I M I I U I 
S 4 0-4-8 

Figure 12: (Case IV) Phase space trajectories in xy plane. 
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2.5 Case V 



Finally let us comment on perhaps the most general form of the solution one could obtain from 
this model, namely that corresponding to turning on all the relevant degrees of freedom. One 
can therefore see that Cases / — IV are actually slices through the full phase space described in 
this section. We will take the following parameterisation for simplicity 



/3 



m 



t4,i 



In this case we will have all three fii non zero which complicates the analysis somewhat, and 
reality again imposes the condition that > 0. Let us initially search for the fixed points around 
z = 0. The primary constraint equation for this becomes 

"-^+f > (2.53) 

Let us initially assume that the denominator is positive definite. Going through the same 
analysis as before yields the usual solution curve = 3x^ provided that we tune /? > and 
a + P > However with reference to the action, we see that this situation leads to both 
W, T diverging and therefore we should be wary of this part of the solution. Returning to the 
constraint equation let us therefore assume that > a + (3 and re-do the analysis. We then 
find that the = 3x^ is perfectly valid, and moreover the parameters W,T are not divergent 
provided that the parameters satisfy a + P — ^ < — (2 + /3). Moreover we also see that /3 is 
bounded from above such that /? < —2/3 - thus severely restricting the form of the variable 
phase space. 

If we search for solutions along the y = slicing things are again complicated. However we 
can simplify things by identifying a = ^, since we can then solve explicitly for x via 

The remaining equation coming from y' = has several solutions. The simplest being = 

give rise to the points 



Xo 



±V2, yo = 0, zq 



xo = ±1, yo = 0, zo = (2.55) 

however the first of these conditions also requires that ^^/^ > Am for the solution to be real. 
The maximal value of ^^/^ is actually given by ^ = which imposes a tight constraint on the 
background parameters which can only be satisfied through substantial fine-tuning. Again more 
general solutions are only available through numeric methods. 
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Figure 14: (Case V) 3-D xyz phase space trajectories for T{(j)) = {(j) / \)°' {(p) = {mcf))^/^ and 
W{<j)) = {(t)/5f. Here, Mp = 1, F = 1, m = 1, A = 1, a = 1, /? = 1, 5 = 1, ^ = 2 and u; = (dust 
case) 






2-1 




-2-' 



Figure 16: (Case V) Evolution of x, y, z versus e-folding number. 
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3 Perturbations and fixed point stability. 



We now need to evaluate the stability of these fixed point solutions. Clearly one may anticipate 
that solutions such as (0, 0, 0) may well be unstable. We must perturb the field equations about 
small values; therefore we need 

x^xo + 6x, y^yo + 5y, z^zq + Sz. (3-1) 

Now the analysis is more complicated than in standard models due to the complexity of the DBI 
action and the general (unknown) phase space dependence of the variables T, W, V. Since 7 is 
independent of any particular parameterisation, we can calculate the general result. 

Using this we can write the perturbation in H' / H. In general we can Taylor expand the function 
W such that we have VK(x* + e*) ~ W{xq) + diWe^ and therefore the general result is true 



1 

1 



+ 



jW \ 3x1 3x1 ^ 



where all terms such as 7, W are evaluated on the classical solution and there is a summation 
over Latin indices. 

The general equations even for the linear perturbation, are shown below for Case V - which 
encompasses all the other solutions in the relevant limit: 

yz^ s f 6y 5z yz^ y^ ( ^^y 



5x = -^(w + /^2)^- + 3--2-J-— (Mm + ^25^2)-^^2■ 



y X J 



.3 



e / f r c- r 1 f ^x 5y\} ^ f ^Sx \ 

Oy = — <IJ-lOf-ll + IJ.20fl2 + m + 1-121 ^ ) J ^ i ~ ^2 ) 



Z^ 



xy 'J yx"^ y x 



y\6y 6x. ^ \ ff^ ^ / H' 



7X \ z \ 6n J \ n x y n 3x'^ \ n 

= '^fiU^^^A.^±^i,\.,A.,i('£\ (3.3) 

2x \ z y X ) Ho \H J 

where we have defined n = a + f3 — p for simplicity and also the following terms 

S^, = 2(a - 2p) Sz {a-2-p) j^ySy ^ 2{a - 2 - p) 6x / ^ j^y^\ 
n z 2n 3x'^ n x \ I2x'^ J 

4(a - 1 - p) 6z {3a-3p-2)fy6y 6x f (3a - 3p - 2)7 V 

op2 = ^— 2 — I I 4(a — I — p) + 



n z n 6x^ nx \ 6x^ 

2{a + 2 + 3p + 2/3) 6z {2a - 2 - Wp + P)'y^y6y 6x {-2a + 2 + lOp - P)j'^y^ 

= ^ a — 2 ^ ^~2 

n z bnx'^ nx bx^ 

5x (2a + 4 + 6p + 4/3) 
x n 
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We will work through an explicit example to illustrate the formalism, namely the Case I solutions. 
Firstly we can calculate the following expression 



1 



1 



3Wxl 



(3.5) 



which will allow us to calculate the perturbed phase space variables. The perturbed dynamic 
expressions then take the following form 



5x' 



by 



5z' 



+ 



z 3x2 



y2 /25y 

2x V y 



5x 

X 



1 



ax 



2x2 



A ^0 JH' 



^(1 

y 



n-f^y 6z 6x nj^y^ . 

' + (3 - n) \ (n - 2 + 



6x2 



z 



X 



6x2 



ax 



3zVi 
2x 



1 + 



(1 



6x2 J fjm^2 
2 

ax 



2^(1 + ^: 

X bx^ 



z 



3x2 



y 



5x , 

X 



6x2 



xy 



[1 



ax 



y^ r^y 



x2 y 



(5xj 

X 



35y(l 



6x^ y 



6x2 



Sxzafii ( 6x 



-[1 +n + 



2Yy 



5z 6y 27 y 



2, ,2 



n. 



y 3x2 (^+4) 



2x 



(52 5x 
- n) 1 n 

Z X 



n72y2 5y 



2 2 
nj y . 

6x2 



- z5 

-no 



(3.6) 



where the notation Hq/Hq implies that we take this function evaluated at the critical points, 
and we have defined n = {a — [3 — 2)/(a — j3) for simplicity. Note that these are the leading 
order solutions only, and that all terms proportional to 5"^ have been neglected. 

The stability of the fixed point solutions is therefore determined by the eigenvalues of the 
resulting perturbation matrix. A lengthy calculation which we will omit here shows that the 
point (0, 0, 0) leads to the eigenvalues 



Ai 



3(6J - 1) 



A2 



3(1 + 0;) 



A, 



3(1 +u;) 



(3.7) 



which indicates that this is never a point of stability for the theory unless the equation of state 
is phantom ie < — 1. In fact this statement will be true for all the various cases we have 
considered in the physical limit, since the dynamical equations of motion all reduce to the exact 
same form in this instance. 



Another relatively simple case to consider is that in Case III. For slices through the (x, y) 
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plane at z = we find the eigenvalues 

A = l{x^ + y^) + l{u:{x^-l)-l) (3.8) 
1 



A± = -^{-Qx^{l + oj) -Qx^ + 2y'^x^ + by'^ ±F{x,y)) 



F{x,y) = ^/l2y^x^io + 96y'^x^io - ASy'^x^uj - 8y^x^ + ASy'^x^ + Idy'^x'^ - 48y^x^ + 36^2x4 + ITy^ , 

If one now slices this through y = we see that we are left with the same situation discussed 
above (as expected), indicative of a phantom equation of state. 

On the other hand, through the y = plane we see that the eigenvalues become 



A = -{l+u;)(l-z'-x'(l-^\ (3.9) 



^(^^ ^(^ 2 2^ Q^^ 
-{l+u:)[l-z -X [l-—^ 

A± = — "^{-Qz^ -x + 2xz^ + x^ -xz^Q + x^ ±F{x,y)) 
2x 

where F is another polynomial in x, 2: and we have defined Q = T/V for simplicity. In the limit 
that z — > 0, we find that these simplify to yield 

A - ^(i+^)(i-x2), 

A± ^ ^(l+w)(x2-l-x2(l±l)). (3.10) 

Note that two of the eigenvalues are therefore degenerate as before, requiring a phantom equation 
of state, however the final eigenvalue has the opposite sign and therefore this fixed point is always 
unstable. 

The remaining fixed points can be analyzed in precisely the same manner, although the 
analysis is somewhat awkward. We will postpone the relevant discussion here and return to it 
in a follow-up publication. 



4 Discussion 



We have initiated an alternate approach to the problem of k-essence, or DBI quintessence |13) . 
using a more generalised form of the DBI action. Since this has more degrees of freedom, 
the resulting analysis is typically complicated, but the phase space structure is far richer. We 
have attempted to make some headway by restricting the phase space volume to various two- 
dimensional slices, and attempting to identify the relevant solution curves upon which the fixed 
points may lie. Our ansatz for each of the unknown functions is also potentially restrictive, 
however we are confident that it represents the leading semi-classical contributions which may 
(or may-not) be derivable from a full string theory embedding of out model. 

What is clear is that the ratio of the (warped) brane tension to the potential is an important 
factor in the dynamics of the theory, where we found T > V in several cases. Moreover the 



30 



additional multiplicative factor W{(p) plays a crucial role, even when it is a constant, since it 
comes into the field equations non-trivially in the expression for H' /H. In the usual DBI analysis, 
W = 1 and the tension is the sole term responsible for the interesting quintessence behaviour. 
In some string compactifications, where the warp factor has no cut-off at small distances, we 
typically find W is constant and greater than unity. However there may be entire classes of 
solution where W < 1, which can lead to novel phase space trajectories. Since our approach has 
been phenomenological, and that there may be additional string backgrounds of interest that 
have yet to be fully explored, we cannot rule out W < 1 - which is vital for obtaining fixed point 
solutions in Case I for example. 

Our numerical results have shown that there is indeed a rich phase space structure present 
due to the increased number of degrees of freedom. We expect many of these to yield highly 
non-trivial stable fixed points in the full analysis, which is beyond the scope of the current note. 
We have classified the nature of as many of the fixed points as is feasible within the current 
analysis. Ultimately we hope that this will lead to a renewed interest in dynamical dark energy 
models driven by a more generalised approach to DS-brane dynamics. 

In light of the recent developments in holographic dark energy [311 E2] and the apparent 
relation to agegraphic [33^ dark energy, we hope that it may be possible to reconstruct the 
various potentials in our generalised model along the lines of |35j . 
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